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A quasi-classical theoretical description of polarization and relaxation of nuclear spins in a quan- 
tum dot with one resident electron is developed for arbitrary mechanisms of electron spin polariza- 
tion. The dependence of the electron-nuclear spin dynamics on the correlation time r c of electron 
spin precession, with frequency Q, in the nuclear hyperfine field is analyzed. It is demonstrated that 
the highest nuclear polarization is achieved for a correlation time close to the period of electron 
spin precession in the nuclear field. For these and larger correlation times, the indirect hyperfine 
field, which acts on nuclear spins, also reaches a maximum. This maximum is of the order of the 
dipole-dipole magnetic field that nuclei create on each other. This value is non-zero even if the 
average electron polarization vanishes. It is shown that the transition from short correlation time 
to fir c > 1 does not affect the general structure of the equation for nuclear spin temperature and 
nuclear polarization in the Knight field, but changes the values of parameters, which now become 
functions of Qt c . For correlation times larger than the precession time of nuclei in the electron 
hyperfine field, it is found that three thermodynamic potentials (x, £, ?) characterize the polarized 
electron-nuclear spin system. The values of these potentials are calculated assuming a sharp tran- 
sition from short to long correlation times, and the relaxation mechanisms of these potentials are 
discussed. The relaxation of the nuclear spin potential is simulated numerically showing that high 
nuclear polarization decreases relaxation rate. 

PACS numbers: 72.25.Rb, 75.75.+a, 76.20.+q, 78.67Hc 



I. INTRODUCTION 

The electron-nuclear spin system (ENSS) of a semi- 
conductor quantum dot (QD) has been under intensive 
investigation in recent years^^. This strong interest 
has been motivated by potential spintronics and quan- 
tum information applications, for which semiconductor 
quantum dots are promising 2,3 . The spin dynamics of 
this system is described by a variety of relaxation times 
which range from nanoseconds to seconds. 

Optical orientation is a commonly used method to cre- 
ate and control^ the ENSS with a high degree of po- 
larization. Nuclear polarization is caused by the Fermi 
hyperfine interaction^ between nuclear spins and photo 
oriented electrons. A simple theoretical description of 
the ENSS behavior is the short correlation time approx- 
imation (SCTA). (The electron correlation time, t c , is 
the characteristic time of the free coherent undisturbed 
electron spin precession in the hyperfine field of the nu- 
clei.) The SCTA is valid if the frequencies of electron 



spin precession (f2) in the local nuclear hyperfine field, 
and nuclear spin precession (u>) in the electron hyperfine 
field are small enough: fir c -C 1, and ujt c C 1. In a 
quantum dot, the electron interacts with a macroscopic 
number, N, of nuclei, i. e., N cx 10 5 and f2 ^> ui. It 
follows that the SCTA can be used if fir c -C 1. 

When f2r c <C 1 holds, frequencies of electron and nu- 
clear spin precession are practically constant during the 
time t c . Small deviations of these frequencies during r c is 
a small perturbations to the spin motion. This deviation 
is the mechanism for the slow transfer of spin polarization 
between electron and nuclei. 

The SCTA is valid in many experimental scenarios. 
However, studying the problem beyond the SCTA stim- 
ulates experimental and theoretical investigations in the 
regime of intermediate correlation time, where f2r c > 1, 
and ujt c 1, and in the regime of long correlation time, 
where lot c 3> 1. These regimes occur at low temperature, 
and under constant wave (CW) light of low intensity or in 
the darkness, respectively. The intermediate regime may 
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be realized also by using circularly polarized light pulses. 
In this case, short periods of light illumination and high 
photoelectron concentration arc alternated with long pe- 
riods of darkness, when the ENSS motion is undisturbed. 

In this paper we discuss the behavior of an ensemble of 
quantum dots, each containing a single resident electron. 
We consider the simplest experimental scenario, where in 
the first step some external action (e. g., circularly polar- 
ized photons) orients the resident electron which polar- 
izes the QD nuclei, and in the second step the ENSS dy- 
namics is determined only by interactions between quan- 
tum dot spins. The electron becomes polarized due to 
the spin exchange with optically oriented photo carriers. 
After switching off the light, a relaxation takes place. 
It is characterized by a long relaxation time T ^> 
The relaxation is a result of the dipole-dipole interaction 
between neighboring nuclei, and of the electron-phonon 
interaction^^. 

In Section [II] we demonstrate (i) that the maximum 
rate of nuclear polarization by optically oriented elec- 
trons is reached for flr c ks 1, (ii) that the nuclear polar- 
ization is a result of the nuclei cooling by spin oriented 
electrons in the Knight field (connected with the time 
averaged electron polarization), and (iii) that the spa- 
tial dependence of the hyperfinc interaction decreases the 
photoinduced nuclear polarization. In the intermediate 
regime flr c mostly influences the nuclear spins relaxation 
times. In Section IIII1 we discuss the difference in the 
ENSS description between the intermediate and long cor- 
relation time approximations. We also present numerical 
results for the dipole-dipole relaxation. These calcula- 
tions show an increasing relaxation time with increasing 
nuclear polarization. All cases arc calculated for a Spher- 
ical Quantum Dot with Infinitely high Barrier (SQDIB), 
which allows us to ignore the exponentially small escape 
of the electron wave function out of the quantum dot. 
We limit the description of the ENSS interaction with 
the QD's environment by introducing the leakage fac- 
tor approximation, and we do not discuss the specifics 
of spin diffusion of the nuclear polarization outside the 
QD&lfiii. 

To model the spin system we use the quasi-classical ap- 
proximation, which is valid for quantum dots with large 
numbers of nuclear spin a 12 i 13 . Finally, in section IIVI we 
summarize our main results. 



nal perturbations of electron field orientation. We do not 
specify the character of the external interaction, but will 
simply suppose that this interaction can partially orient 
the electron spin. The exchange scattering of an opti- 
cally polarized carrier is an example of such an external 
interaction. 



A. Electron spin precession in the nuclear 
hyperfine field 

In GaAs-like semiconductors the electron and nuclear 
spins are coupled by the Fermi hyperfinc interaction, 



(1) 



where \xb is the Bohr magneton, s and r are the spin and 
position of the electron, fj, n , I n and R n are the magnetic 
moment, spin and position of the n-th nucleus, respec- 
tively. The sum in Eq. {TJ) runs over all the nuclei inside 
the QD. The hyperfine energy has a maximum when s 
and I n are parallel, and a minimum when they are anti- 
parallel. In the following, and for simplicity, we will sup- 
pose that all nuclei have the same spins and magnetic 
moments I n = I and \x n = fj,i. 

A QD contains a macroscopic number, N, of nuclear 
spins, N oc 10 5 ^> 1. Therefore, the frequency of electron 
spin precession in the nuclear hyperfinc field, 



ft 



(2) 



is distributed in a wide region from zero to f2 m ax given 
by 



(3) 



It is usual to separate ft in two parts, average and fluc- 
tuation: ft = (ft) + Aft. In the following, average and 
fluctuation refer to the time evolution and time fluctua- 
tion of the nuclear spins. Because (ft) and Aft are not 
correlated, (O 2 ) = (O) 2 + ((AO) 2 ). 

ft is many orders of magnitude larger than the fre- 
quency of nuclear spin precession in the electron hyper- 
fine field, u) n s, i. e., ft > flfi uc , where 



II. POLARIZATION OF THE 
ELECTRON-NUCLEAR SPIN SYSTEM IN A 
QUANTUM DOT 

In a semiconductor quantum dot with one resident 
electron the hyperfine interaction creates a localized 
electron-nuclear spin polaro n 14 ' 15 i 16 . In this section we 
discuss this system's behavior in the limit of short and 
intermediate correlation time, ujt c <C 1. In these regimes 
the frequency of the nuclear spin precession in the elec- 
tron hyperfinc field is lower than the frequency of exter- 
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Here, flfi uc = \J (ft 2 ) \au=o is the characteristic value of 
the fluctuations of the electron spin precession frequency, 
|/|| = \Jl(I + 1), and = \J s(s + 1), arc the modulus 
of nuclear and electron spin, respectively, 
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and ip(R n ) is the electron wave- function on the n-th 
nucleus^. 

^max does not depend on the quantum dot's volume 
because the electron wave function is normalized in this 
volume; it is determined solely by the chemical compo- 
sition of the QD. For example, for a GaAs QD one can 
estimate^ f2 max ~ 10 _11 s _1 . Values of f2fl uc oc ri max /v^V 
and (lu) oc fl max /N depend on the QD volume. For a 
quantum dot with N oc 10 5 nuclei fla uc oc 3.10 8 s _1 and 



(lu) oc 10 6 



Characteristic frequencies and times arc 



collected in Fig. Q] 

Under sample illumination, the photocarriers and pho- 
tons scattering on a quantum dot, carrier capture and 
photon absorption in QD are the main mechanisms of 
the ENSS interaction with the environment. As we will 
show in section III BL the nuclear spin polarization and 
relaxation are a consequence of this perturbation of the 
electron and nuclear spin precession in the hyperfinc field. 
To achieve the most effective nuclear polarization, the fre- 
quency of these collisions, lu co u = t~ x , should be about 
fl 3> lu h . Between collisions nuclear spins change their 
directions by a very small angle of about t c (lu)/2, which 
is much less than 1. In the zero-th order approximation, 
the spin, s(t), of the resident electron precesses in the 
constant (frozen) nuclear field: 



s(t) 



sn(h) + (s(t ) - s n (t )) cos(fl ■ St) 

+ [en x s(t )]sin(fJ • St), (6) 



where 



sn(to) = (s(t ) ■ e n ) ■ e n , 



(7) 



efj = £~2/r2, s(to) is the initial spin, which is determined 
by collisions, St = t — to, and to is the time when the 
external action polarized the resident electron. 

For the continuous-wave (CW) photoexcitation, the 
probability W(St) of electron free (undisturbed) preces- 
sion during time St decreases exponentially with free pre- 
cession time: W(St) = exp(— St)/r c . The spin of localized 
electrons (averaged over initial polarization s(to) and St) 
can be written as: 



so 



[fl X S ] T c + (fl ■ S )flT^ 



1 + (fir c y 



(8) 



where Sq is the average value of the initial spin. For an 
excitation with a periodic train of short pulses W(t) ss 
S(t — t — t c ), and s is given by Eq. ©, by substitut- 
ing cos(fl(t — to)) — > sm(flr c ) I (flr c ) and sin(f2(t — to)) ~> 
(1 — cos(flT c )) / (flr c ) . For our goal, the difference between 
the CW and the pulsed excitation regimes is only quan- 
titative. In this paper we analyze (and present model 
calculations) only for the CW excitation regime. 

In the short correlation time approximation s w So, 
whereas in the opposite limit (flr c > 1), the mean value 
of the electron spin depends on the angular distribution 
of frequencies, fl. For a random distribution (unpolar- 
ized nuclear system) (fl) =0, (fl ■ sq) ■ fl = il 2 So/3, 
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FIG. 1: Time scale for different regimes of hyperfine inter- 
action in quantum dot. (1) Short correlation time regime 
(t c < £l~ 1 ). During r c the electron and nuclear spins ro- 
tate around a small angle. (2) Intermediate correlation time 
regime < r c < (lu)^ 1 ). During r c the electron spin 

rotates around a large angle with constant angular velocity. 
(3) Long correlation time regime (r c > (a;) -1 ). During r c 
the electron and nuclear spins change their direction. For 
(<^dd} _1 > t c > (cu)^ 1 the total nuclear spin, It, is conserved, 
but Q changes direction as a result of the nuclear spin preces- 
sion in the non-uniform electron hyperfine field. (uJdd) is the 
frequency of nuclear spin precession in the local magnetic field 
of the neighboring nuclei. For (ujdd)" 1 < t c , the dipole-dipole 
interaction between spins of neighboring nuclei changes 7e- 



and s ss So/3. For a polarized nuclear spin system 
(fl) ^> flfiuc, and the mean value of the electron spin 
depends on the direction of (fl). In the absence of exter- 
nal magnetic fields the ENSS has only one distinguished 
direction that is determined by the photo-electron polar- 
ization, Sq. For (fl) along Sq (fl ■ sq) ~ flso, and s w Sq. 



B. Nuclear spin precession in the electron 
hyperfine field 

Let us now describe the mechanism of the nuclear po- 
larization by the resident electron. We consider a slow 
nuclear spin precession that obeys the equation 

^ = L0 n [s(t) X I n ] 

at 



UU n \S X J„ 



ds 
~dt 



x J„ 



St 



S X 



din 

dt 



St. 



(9) 



The average in this equation is done on the time region 
At, such that fl <C At -c lu. The first term in the 
right-hand side of Eq. © gives the regular part of the 
nuclear spin precession in the mean electron hyperfine 
field, Bpc n = —LUn&, known as the Knight ficld^i. The 
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second term describes the dynamical polarization of the 
nuclei. For small nuclear polarization, (/) <C /, and: 



din | _ / [[n x s ] x I n ] ' 
dt ldp ~ W " Tc \ l + (fir c )2 fl - (^)|| s ||2 Tle 



where 



ri e (n) 



211 -w I 



(fir c ) 



3(1 + (fir c ) 2 ) 



(10) 



(11) 



is the characteristic time of the longitudinal nuclear spin 
relaxation, and (uj 2 ) = ^2 n u n /N. The third terrain in 
Eq. ®, 
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din 

dt 



^]w m s(t) x / [s(ti) x I n ]dh 



(12) 



represents two processes: relaxation of nuclear polariza- 
tion (averaged over initial electron spin direction), 



dl, 



" 2 Jn , 2+(Or e )% 

5i ke '~ (^)T le (fi,r c ) ( + 2 Udj 



and nuclear spin precession, dl n /dt\i nc t = [t7„ x I n ], in 
the indirect hyperfine field given by 



4 [n x i n ] r c 

^)T le {n,T C ) 



(14) 



(For a derivation of Eqs. (TO" ]) . ([T3" j) and (fHj) see the 
appendix.) Both of these processes, i. e., relaxation of 
nuclear polarization and nuclear spin precession, are de- 
termined by fluctuations of the electron hyperfine field. 
They are proportional to ||s|| 2 = 3/4 because T\ e oc 
| t s||~ 2 . The relaxation time versus flr c is presented in 
Fig. [21 In the limit of short correlation time, all com- 
ponents of nuclear polarization relax at the same rate*, 



T^O,^)" 1 =2|| S || 2 T C 



2 )/3. 



For intermediate correlation times, the relaxation time 
for the component of /, In, longitudinal to ft, increases 
as both f2 and r c increase, yielding^ Ti e (fi, r c ) = (1 + 
(flT c ) 2 )T le (0,T c ) (solid curves 1,3 in Fig.© On the other 
hand, the relaxation rate of the polarization compo- 
nent, I±, transverse to ft, behaves different depending on 
whether Q is increase or r c is increased. For t c increasing, 
Ti e ,J_ decreases monotonically (curve 2). As il increases, 
T\ c ,x increases saturating at Ti e ^(oo,T c ) = 2Ti e (0,T c ) 
(curve 4)2S. 

The nuclear polarization, its relaxation rate, and the 
indirect hyperfine field, are all proportional to uj 2 t c oc 
HV^-Rn)!! 4 , and all have a strong spatial dependence. 
They have a maximum in the center of the QD, and they 
decrease towards the border. In the limit of frozen nu- 
clear field, i. e., short and intermediate correlation times, 
u> n T c in Eqs. (fTO)) . JT3]) and (fTl)) is much smaller than 



FIG. 2: Nuclear spin relaxation time on the resident elec- 
tron vs. correlation time r c (curves 1,2), and vs. frequency 
Q (curves 3,4) of the electron spin precession in the nuclear 
hyperfine field. 

Solid curves 1,3 show the relaxation time for the nuclear polar- 
ization component along ft. Dashed curves 2,3 show the relax- 
ation time of the nuclear polarization component transversal 
to f2. Dependences of relaxation time on r c (1,2) were cal- 
culated for a constant value of fi. They are normalized to 
the value of Ti. c at fir c = 1. Dependences of relaxation time 
on 51 (3.4) were calculated for a constant value of r c , and are 
normalized to the value of Ti, c at = The nuclear dynamic 
polarization is ineffective for f2r c 3> 1 because in this region 
the value of the relaxation time for the nuclear spin longitu- 
dinal component increases fast, and the leakage factor in the 
Eq. (|16[) decreases. The difference between times Ti >e for lon- 
gitudinal and transverse component of nuclear polarization 
in the region Qt c 2> 1 increases e in Eq. (|16|l . but decreases 
nuclear polarization. 



1. In Eq. (|14p . for an intermediate time, this parameter 
may be rewritten as w n /fi, which is much less than 1, 
implying that characteristic rates for nuclear spin polar- 
ization, relaxation, and indirect hyperfine interaction are 
many times less then tu n - These slow processes are rele- 
vant only if they change the system's behavior, such as 
its nuclear polarization, Eq. (|10|) . or its relaxation rate, 
Eq. CE2). 

By comparing the indirect field and the Knight field 
one can see that the former is significant only for sq < 
(w„/fi) • (flr c ) 2 1 (1 + (Qt c ) 2 ) < N- 1 / 2 < 1. In the follow- 
ing context we suppose that the electron polarization is 
high enough, which allows us to ignore the indirect hyper- 
fine interaction when it comes together with the Knight 
field. (The indirect hyperfine interaction between nuclei 
plays an important role in the relaxation of the electron 
polarization transversal to a strong external magnetic 
fieJdiLliia). 

A sketch showing the main mechanisms of the hyper- 
fine interaction's influence on the ENSS's behavior under 
sample illumination is presented in Fig. [3l This figure 
shows three precessions (i) in the mean Knight field, cre- 
ated on the nuclei by a mean electron polarization, (ii) in 
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clear polarization, T± e 
reached when r. 



(min) 



10" 
= 10 



-11 , 



This relaxation time is 
Since the short corre- 



External mechanism of the resident 
electron polarization 




FIG. 3: Different mechanisms for the hyperfine interaction 
between the resident electron and the nuclei. 
Nuclear spins precess in the mean Knight field, and are po- 
larized by the oriented resident electron. The rates of these 
effects are proportional to the average electron spin. Polar- 
ized nuclei create a mean Owerhauser field on the electron, 
and a mean indirect Weiss field on each other. Their polar- 
ization also relaxes via interaction with the resident electron. 
The rates of these three mechanisms are proportional to the 
value of nuclear polarization. 



the mean Ovcrhauser field, created on the electron by the 
mean nuclear spin, and (iii) in the mean indirect Weiss 
field, created on the nuclei by the mean nuclear polariza- 
tion. It also contains two dissipation processes: nuclear 
polarization by the oriented electron through the fluctua- 
tion of nuclear spins, and nuclear polarization relaxation 
on fluctuations of the electron hyperfine field. 

When transitioning from short to intermediate corre- 
lation times, what changes most is the dependence of 
the Ti c time on Qt c . For f2r c <C 1, the time T± e is 
proportional to (cj 2 ) -1 , and also to TV 2 . T\ e is a de- 
creasing function of r c , and it reaches its minimum for 
Qt c = 1. On the other hand, in the limit Qt c 3> 1, 
T le (tt) pa T c (3fJ 2 )/(2||s|| 2 (cj 2 )) is an increasing function 
of t c . For a depolarized nuclear system, i. e., when 
O pa ttfiuc, the previous equation gives T\ e cx t c N, 
whereas for a polarized nuclear system T\ e cx t c N 2 , as 
we saw before. 

For a typical GaAs QD (with N = 10 5 , (w) = 10 6 s -1 , 
f!g uc w 3 ■ 10 8 s _1 ), the shortest relaxation time for de- 
polarized nuclei, calculated from Eq. (|11[) and condi- 
tion dTi e /d,T c = 0, 

(min) 



3 • 10" 



is Ti e ~ 10 3 s. It is reached for 
On the other hand, for 100% nu- 



lation time approximation is valid for r c <C Tq 111111 ' . this 



(min) 

c j 

region decreases by a factor of VN for a highly polarized 
nuclear system. 

The relaxation rate and precession frequency depend 
on the nuclear position in the quantum dot. The balance 
between nuclear dynamical polarization, Eq. (| 10[) . and 
relaxation, Eq. (fTB")) . gives: 



|/1| 2 S0 

INI 2 ' 



(15) 



Since the right-hand side of this equation does not con- 
tain n or t c , the average nuclear spin has the same value 
for all nuclei, i. e., (I n ) = (I). The previous statement is 
valid both for short and intermediate correlation times. 

Equation (p~5|) is correct if we take into account only 
one mechanism of nuclear spin relaxation, namely the hy- 
perfine interaction with fluctuations of electron polariza- 
tion. Additional channels of relaxation decrease the nu- 
clear polarization. Often this decrease can be described 
by introducing in Eq. (|15[) the phenomenological leakage 
factori / = T u /(T le + T u ) < 1. (Here T u is the relax- 
ation time for additional relaxation channels.) When Tu 
is finite, the leakage factor, /, and the average polariza- 
tion, (J), are monotonously decreasing functions of T le , 
and have a maximum for flr c pa 1. 



C. Dipole-dipole interaction between nuclear spins. 
Nuclear spin temperature 

The main additional channel of nuclear polarization 
relaxation is determined by the dipole-dipole interac- 
tion between neighboring nuclear spins. This interaction 
transfers nuclear angular momentum to the crystal lat- 
tice with a characteristic time T2 ~ 10~ 4 s, which is much 
smaller than^ T\ e . 

As a result, the steady state value of the quasi equilib- 
rium nuclear polarization is TijT\ e times less than that 
predicted by Eq. (|T5)) . and it has to be reached at a 
time T2. Nevertheless, it is well knowni^ from many 
experiments and theoretical calculations, that the op- 
tically induced nuclear polarization cannot usually be 
ignored. This polarization is due to a decrease of the 
nuclear spin temperature, O. The effect results from 
the balance of two energetic flows: cooling of nuclear 
spins by oriented electrons in an external magnetic field 
B. J coo i cx — (B ■ <9I/di|d P ), and heating of nuclear 
spins by random fluctuations of electron polarization: 
Jheat = PC Here (i = (/ss©) -1 is the inverse spin 
temperature, and C = dE/d(3 is the heat capacity of 
the nuclear spin system. In the short correlation time 
approximation, (f2r c <^ 1), and in a spatially uniform 
external magnetic field^ 



/? = / 



4/ (B-a ) 

m s 2 + eBi ' 



(16) 
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where B\ is the characteristic value of the random local 
field squared, and e is a number of the order of one. j3 
has the opposite sign for parallel and anti-parallel orien- 
tations of B and Sq. (It is not surprising that the spin 
temperature can be negative, since the energy of the spin 
system is limited both from above and below^^). (3 is 
positive if Sq is parallel to B, in which case spin align- 
ment decreases the nuclear spin energy. 

In a general case, the field Bl is a result of 
dipole-dipole and indirect hyperfine interactions between 
nuclei^. For the special case of Eq. p4[) . the indirect 
interaction between two nuclei, n and m, depends on n 
and m only through their product, u! n u! m , and has no 
influence on the local field part of Eq. (fT6|) : eB\ rts 3B 2 d . 
Here B dd is the square of the average dipole-dipole part 
of the local field. 

Eq. (|16p is valid for high nuclear spin temperature, i. e., 
for fi\xi\J B 2 + eB 2 L <C 1, and was derived for a spatially 
uniform B and T\ e . An average Knight field should be 
included in the regular external field, and the total mag- 
netic field, B + Bfc,m depends on the nuclear position 
inside the QD. 

The approximation that ignores the spatial depen- 
dence of the hyperfine interaction inside the dot is 
know n 21 ! 22 as the "box model". For the "box model" 
the nuclear spin temperature and nuclear polarization 
have the same value for all nuclei. In the real situation 
of a spatially inhomogeneous hyperfine interaction, the 
nuclear spin temperature and polarization cannot both 
be constant because (J n ) oc (3Bk, n - 

The dipole-dipole interaction between nuclear spins 
produces an energy flow from the region with high spin 
temperature to the region with low spin temperature. It 
is commonly assumed that the nuclear spin diffusion in- 
side the area of electron localization is suppressed by the 
strong gradient of the Knight field. (If the difference in 
hyperfine splitting of the nearest nuclear spin levels is 
larger than their dipole-dipole broadening, the flip-flop 
process between nearest nuclear spins is suppressed by 
the energy conservation law.) But one can show that 
the typical difference in the splitting of the nearest nu- 
clear spin levels for a spherical QD is about h(ui)so/Nji ft 
M7 max so(4// 3 iV 4 ) 1/3 , where N R ft (iV/4) 1 / 3 is the num- 
ber of nuclei along the QD radius. For a GaAs QD 
with N = 10 5 this difference is about or less than 
hl0 4 s~ 1 ft hoJdd, and, therefore, there are no reasons 
for a strong suppression of the spin diffusion. 

In the limit of efficient spin diffusion, the nuclear spin 
temperature should have the same value for all nuclei, 
and the equation for (3 contains averaged values, i. e., 

p = -± f (17) 

Here (uj m ) = £u;™/iV and = A B \I^ 2 - 

n 

In the limit of short correlation time, the differ- 
ence between the result of Eq. (JT7J) and the one ob- 
tained with the "the box model" is only numerical. 



The dimensionless saturation value of the spin tempera- 
ture, (3 = h(uj}(3, is then {lu) {u 3 ) / (lu 4 ) times less than 
that predicted by Eq. (fT6|) . It reaches saturation if 
Sq 2> e (lj dd ) (lu 2 ) / (lj 4 ). (For a spherical quantum dot 
with infinite barrier (SCDIB), we have: (lu) 2 / (lu 2 ) ft 
0.36, (lu)(lu 3 )/(lu 4 ) ft 0.22, (lu 2 ) 2 / (lu 4 ) ft 0.17, and 
£ ( W 2 d )(c 2 )/(^ 4 )«5-10- 3 ). 

In equilibrium, the mean value of the nuclear spin in 
this field is 

(I n > = -/Ji(fe n ) So . (18) 

When there is a single global nuclear spin temperature for 
all points in the QD, the maximum average nuclear po- 
larization within the SCDIB model is 4 to 5 times smaller 
than the one predicted by the "the box model" . 

In section Hi Bl we saw that T le as function of r c has a 
minimum for t c = il^ 1 . For t c > f2 , increasing the cor- 
relation time and the frequency 1~2 increases T± e . The nu- 
clear polarization increases H, implying that if one wants 
to optimize the nuclear polarization by varying the cor- 
relation time, one has to start from a unpolarized nuclear 
system with fiflucT'e C 1. At the same time, r c should 
be long enough to achieve the condition f2r c = 1, and to 
maximize the leakage factor in the final polarized state. 
Then, the best regime for generating high nuclear polar- 
ization is on the border between short and intermediate 
correlation times. 

For flr c 3> 1, the anisotropy of the nuclear spin re- 
laxation on the electron also renormalizes the param- 
eter £ s=! 3 + (f2r c ) 2 . We can neglect this effect for 
Qt c < 1. In Fig. [U results of the nuclear spin polar- 
ization calculations for the "box model" (curve 1) and 




o.o 0.2 0.4 

S o 

FIG. 4: Dependence of nuclear polarization on the mean spin 
of the resident electron. Here Ti e (0) = Tu- Curve 1 was 
calculated using the "box model" with fi m axT c <g 1, curve 2 
using the SCIB model with Q max rc <C 1, and curve 3 using 
the SCIB model with f2 ma xT c = 10. 
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SQDIB model (curves 2,3) are presented. All calcula- 
tions are done for T le (0) = T\;. Curves 1 and 2 are 
calculated in the short correlation time approximation, 
curve 3 for £l m & x T c = 10. One can see that the nuclear 
polarization for the SQDIB model is about 5 times less 
than for the "box model" . For f2 max T c = 10, increasing 
the time T\ e by increasing the nuclear polarization addi- 
tionally decreases the leakage factor and slows down an 
increase of the nuclear polarization. 



III. ELECTRON-NUCLEAR SPIN SYSTEM IN 
THE LIMIT OF LONG CORRELATION TIME (IN 
DARKNESS) 

In this section wc consider the relaxation of an isolated 
ENSS in the limit of long correlation time. The nuclear 
polarization by spin-oriented carriers is extremely ineffec- 
tive. We will consider only spin relaxation of an isolated 
ENSS in the darkness. In section IIII Al we discuss the 
main thermodynamic potentials that characterizes the 
isolated spin system. We connect these potentials to the 
relaxation of ENSS parameters under illumination and 
discuss the relaxation mechanism. In section IIIIBI wc 
present the results of a numerical simulation of the QD 
spin relaxation due to the dipolc-dipole interaction that 
transfers nuclear spin into the crystalline lattice angu- 
lar momentum. This process is controlled only by the 
state of the QD ENSS, and is independent of the QD's 
environment. 



A. Conservation laws and thermodynamic 
potentials of the ENSS 

The r c of the ENSS largely increases in the darkness. 
At liquid helium temperatures the characteristic time of 
spin relaxation for electrons on phonons is about seconds, 
whereas for nuclei it ranges from days to yearsi^i^. The 
direct transfer of spin angular momentum to the crystal 
by the dipole-dipole interaction between nuclear spins is 
the main mechanism for the relaxation of the spin polar- 
ization. This interaction is also responsible for the en- 
ergy diffusion from the quantum dot to the neighboring 
nuclei2^25. 

In the zero-th order approximation, we keep only the 
hyperfinc interaction and switch off all other interactions, 
which makes it easier to determine some integrals of mo- 
tion. The Fermi interaction conserves energy, E = HQsq, 
and total spin, F = Ir, + s. The total spin of the nu- 
clei, Je, is many orders of magnitude larger than that 
of the electron. Therefore, with high precision 1-^ is also 
conserved. Moreover, as a result of the adiabatic approx- 
imation, Q 3> Lu n , and the electron spin projection, sq, 
along fl is also conserved and must be quantized, i. e., 
sq = ±1/2. The conservation of follows from the con- 
servation of energy and sq. 

Under illumination every nucleus is acted on by an 



average Knight field, whose value and direction is deter- 
mined by the average electron spin, whereas in the dark- 
ness the ensemble of quantum dots decomposes in two 
sub-ensembles (nuclear spin polarons) with defined value 
of electron spin projection, sq, on the nuclear hyperfine 
field, i. e., sq — ±1/2, and energies ±M7/2. The proba- 
bility to find a quantum dot in one of these sub-cnscmbles 
is given by the conservation of energy, frequency, and to- 
tal nuclear spin. 

We will assume that the transition from illumination 
to darkness is sharp, and that in the initial state the total 
nuclear polarization is larger than its fluctuations, i. e., 
{I) 1/ \N. Considering the integrals of motion of the 
system, we can write the ENSS probability distribution 
$, as a function of three thermodynamic potentials: an 
electron spin potential, ^, an inverse nuclear spin tem- 
perature, Xsn-i f° r each sub-ensemble, and a nuclear spin 
potential, £(x TO ), 



exp{(^s^) - (xsgfiilsn) + (g • Je)} 
(47r)*(exp{?/2} + exp{-^/2}) 



(19) 



At the last moment of illumination, the total nuclear 
spin mean value, (Is), and the electron spin, Sq, are di- 
rected along the Knight field, and the nuclear inverse 
spin temperature under illumination is (3. Therefore, af- 
ter switching off the light, the nuclear inverse spin tem- 
perature in the darkness is 



Xsn « 0— ■ 
SQ 



(20) 



Moreover, nuclear and electron spin potentials are given 
by 



3</y 



and 



? = hi 



1^1 



1 + 25Q 

l-2s ' 



(21) 



(22) 



In Eq. ([20]) wc took into account that for a cooled QD 
nuclear spin system (||/ s || 2 ) « ||I|| 2 A + ||(/s) 2 ||. The 
first part of the right-hand side in this equation describes 
the fluctuation of the total nuclear spin, and the second 
part the square of the mean value of the total nuclear spin 
in the electron hyperfinc field. Equations (|20|) . (HH), and 
(|2"2"|) completely determine the initial state of the system. 
From Eq. ([20]) one can see that X+1/2 = — X-i/2- In the 
following we will consider only the sub-ensemble with a 
positive spin temperature, and omit sub-index. 

In the darkness, the relaxation of the potentials £, Xj 
and is due to the dipole-dipole interaction between 
neighboring nuclei, and to phonon scattering. Each re- 
laxation potential has its own relaxation time, for £, 
for Xi T x ; and for <;, T s . The dipole-dipole interaction 
does not conserve total spin, 1^. As a result, the angu- 
lar distribution of Is tends to the isotropic distribution, 
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Nuclear spin 
Temperature 




Hyperfine 
Energy 

h\\ y. n u) n I n \\sn 



phonons 
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Neighboring 
Nuclei to QD 



FIG. 5: Relations between the ENSS's description under illu- 
mination and in the darkness. Under illumination the ENSS 
state is described by the average electron spin, so, and the 
inverse nuclear spin temperature, j3. The Knight field of the 
electron creates the mean polarization of the cooled nuclei. 
After a sharp transition to darkness, the ensemble of quan- 
tum dots splits in two sub-ensembles with two different elec- 
tron spin projections, sn = ±1/2, on the nuclear hyperfine 
field. The ENSS is characterized by an electron potential, 
?, by a nuclear potential, £, and by a nuclear inverse spin 
temperature, x«n • The relaxation of these thermodynamic 
potentials is due to the dipole-dipole interaction between nu- 
clear spins and electron phonon interaction. This relaxation 
transfers nuclear spin to the crystalline lattice (on a time T$), 
energy diffusion from the quantum dot to the environment 
(on a time T x ), and energy from the scattering of phonons 
(on a time T q ). 



and £ tends to zero. But the dipole-dipole interaction 
conserves total energy, and the modulus of Js fluctuates 
around its average value (Is)(x) = X^( w )-^ll^l| 2 /6- In 
section UlI Bl wc demonstrate that the rate at which £ re- 
laxes depends on the value of Is, or, in other words, that 
is a function of \. 

The relaxation of Js , and of the frequency, (f2) , is a 
result of the nuclear spin energy flow out of the quantum 
dots provided by spin-spin interactions between nearest- 
neighboring nuclei. As a result of this process, x tends to 
0, and the frequency of electron spin precession relaxes 
from its initial value to an asymptotic fluctuation value 
given by f2fl uc = \/||/|| 2 iV(ti; 2 ), which is non zero for a 



finite system. The rate of this relaxation process depends 
on the value and sign of the nuclear spin temperature 
around the QD. 

The state with x ~ an d ? ^ can be identified as 
a fluctuating nuclear spin polaroid. In the fluctuation 
polaron state the electron "remembers" its spin direction 
after total nuclear spin relaxation takes place. To put 
it another way, the potential ? ^ describes not elec- 
tron polarization but correlation between electron spin 
and nuclear field. The relaxation of <j is caused by the 
flipping of the electron spin. The energy of this transi- 
tion, Ml, cannot be taken from the nuclear spin system 
as il ^> ujl, but is due to phonon-scattering instead§»L£. 
A sketch of connections between ENSS parameters and 
their relaxations is presented in Fig. O 

In short, the relaxation of 5 and x 1S due to the open 
character of the QD ENSS, i. e., it depends on the envi- 
ronment. We will not discuss it later. 



B. The dipole-dipole relaxation of an isolated 
quantum dot 

The dipole-dipole relaxation of the nuclear spin sys- 
tem is commonly^ characterized by a time T 2 
For GaAsi T 2 « l(T 4 s. The estimations in Refs.[l||2£ 
demonstrate that, for a QD composed of an electron and 
polarized nuclei (nuclear spin polaron), the relaxation 
time of spin polarization should increase as 



(23) 



where = (1^) /N 2 \\I\\ 2 is determined by the nuclear 
spin temperature, x- Iq is the projection of the total nu- 
clear spin on the direction of the nuclear hyperfine field. 

We now compare the estimation of Eq. (|23|) with the 
results of our numerical simulations of the ENSS spin dy- 
namics. For this we calculate the nuclear spin correlator, 



G{t,pa) 



/I E (t')I s (f' + f)-dt 
/I s (tf).J E (tf)dt ■ 



(24) 



Our numerical model considers a spherical quantum dot 
containing AT mo( j nuclei located on a cubic crystal lat- 
tice. To make the calculation feasible we take N mo< i to 
be of the order of hundreds, and, therefore, much less 
than the number of nuclei in a real quantum dot. How- 
ever, iVmod 3> 1, and the adiabatic approximation may 
be used to model the system. Therefore, the electron 
spin has a constant projection, sr>, on the total nuclear 
field given by sq = ±1/2, whereas the nuclear spin pro- 
cesses around the hyperfine electron field (directed along 
ft = J2 n u nln) with frequency w„/2. 

As the frequency unit we take the mean frequency of 
nuclear spin precession, (u>)/2 = 1. As the length unit 
we take the distance between nearest nuclei, r = 1. We 
use the following equation to calculate the precession of 
the nuclear spin at site n in the magnetic field created 
by neighboring spins, 
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din, 

dt 



dd 



rn^n 



[{\l m - e(el m )) x I n ] 



rf 3 



(25) 



Here e is the unit vector joining sites n with spin to, 
d m „ is the distance between these sites, 7 = u)£>/(u>) « 
5.7 x 10 -2 , and ujd = ^ n /(^ r o) is the characteristic fre- 
quency of the nuclear spin precession in the field of near- 
est neighbor dipoles. 

The sign of the electron spin affects only the direction 
of the nuclear hyperfine precession, but has no influence 
on the dipole relaxation. For this reason, we simulated 
only the sub-ensemble with positive spin temperature. 
The initial distribution of nuclear spins is determined by 
the inverse temperature, Xj with natural units H(lu). This 
initial spin distribution was generated by a random pro- 
cess using the Boltzmann distribution function. In order 
to decrease the effects of crystal magnetic anisotropy in 
the initial stated, f2 was directed along the [111] axis. 

The calculated dependence of pa on the inverse spin 
temperature, x, is presented in Fig. [51 The numerical re- 
sults are in a good agreement with the simple theoretical 
equation, 



(pa) * JW+WJWW))- 



(26) 



Here (L) is the Langevin function, averaged on the QD 
volume, and is given by 



(L(x)) 



sin (r) 



r 2 dr, 



(27) 



where x = x w 0j and the index indicates the nucleus at 
the center of the QD. The second term under the square 
root in Eq. ([26]) describes the nuclear spin fluctuation, 
which is important in the limit of extremely high tem- 
perature, when x^( w ) ^ V v^V- 

To control the numerical precision, we checked the 
conservation law for the energy and the total nuclear 
spin of the system without dipole-dipole interaction (see 
Fig. EJa)). After introducing the dipole-dipole interac- 
tion (Fig. [7Kb)) the total energy is still conserved (curve 
3) , but now all three components of the total spin have a 
random behavior (sec, for example, curve 1 for Js iX (<)). 
The total nuclear spin, fluctuates around its 

mean nonzero value that is related to the total hyperfine 
energy, (to)(I) w (f2) oc E. The fluctuations of hyper- 
fine energy are negligibly small because they reflect very 
weak transfers of energy from hyperfine to dipole-dipole 
reservoirs and vice versa. 

In Fig. M is the calculated correlator G(t), Eq. ([24]) . 
for various nuclear spin temperatures, as indicated. One 
can see that an increase of the nuclear polarization, from 
curve 4 with po = 0.04, to curve 1 with po = 0.152, de- 
creases the spin relaxation rate. On a short time scale, 
the time dependence of the spin correlator can be ap- 




0.1 0.2 0.3 0.4 0.5 

FIG. 6: Nuclear spin polarization as a function of the inverse 
nuclear spin temperature for a spherical quantum dot with 
infinite wall. The solid curve is the averaged Langevin dis- 
tribution, Eq. (|27|l . for infinite numbers of nuclei, the dashed 
curve, the calculation for the SQDIB model, Eq. (|26|l . with 
TV = 489. Circles are the Monte Carlo simulation, also with 
TV = 489. 
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FIG. 7: Conservation of the quantum dot's parameters, (a) 
With the dipole-dipole interaction switched off, the spin com- 
ponents Is,x, Is, y , and Is,z, given by curves 1, 2, and 3, 
respectively, the spin modulus, Is, given by curve 4, and the 
energy of the system, E, given by curve 5, are all conserved, 
(b) When the dipole-dipole interaction is switched on, the 
spin components are no longer conserved, and change their 
values chaotically instead (curves 1, 2, z component is now 
shown). But the spin modulus, curve 3, is still conserved 
within numerical error. 



proximated by a Gaussian distribution: 



G(t, (7),iV)«exp. 



THp,N) 



(28) 



This equation allows us to quantitatively compare the 
result of our numerical experiment with the estimation of 
Eq. ([2"5]) . Figure[5]shows our numerically calculated value 
for T^(/3,N) for 4 different spherical quantum dots, with 
N equal to 251, 485, 895 and 1365 nuclear spins. There 
T^(/3, N)(ui)/N vs. p^fl) is presented; the polarization 
relaxation rate decreases fast with increasing nuclear po- 
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t 1 1 1 1 1 1 1 1 1 1 1 r 




FIG. 8: Spin correlator, G, vs. time for different nuclear tem- 
peratures, as indicated, on a quantum dot with N = 1365 
nuclei. Solid lines are for the numerical simulation, whereas 
dashed lines are a Gauss approximation with fitting param- 
eter T ( following Eq. (28J). (1) /3h{u) « 3.5, pa = 0.410, 
T $ (oj) = 820 ± 15, (2) f3h{uj) w 0.18(p n « 0.051, T $ (oj) = 
160 ±15, (3) f3h{u) « 0.078, pn = 0.036, T<<w) = 90 ±15, (4) 
j3h{w) « 0.05 pn = 0.027, T € (w) = 60 ± 15. 
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that Ti e oc t c . The maximal rate of nuclear polarization 
by the QD's electron is reached for an intermediate value 
of correlation time, Q,t c sa 1. In this case the leakage 
factor reaches its maximum. Nuclear polarization in the 
long correlation time regime is not efficient, which fol- 
lows directly from the general equations in Ref. [28l . that 
connect nuclear polarization and relaxation rates with 
electron spin correlator. 

(2) The nuclear spin diffusion inside a QD increases 
the average nuclear spin temperature and decreases nu- 
clear polarization. The diffusion is a result of the spatial 
dependence of the spin relaxation rate and of the Knight 
field. Nuclear polarization diffuses from the QD's cen- 
ter to its periphery. In the vicinity of the barrier, the 
Knight field is comparable with the local dipole field, and 
strong dipolc-dipolc relaxation destroys nuclear polariza- 
tion. This effect decreases the mean value of the nuclear 
polarization by a factor of more than 4. One should take 
this effect into account when experimentally describing 
the value of the nuclear polarization. 

(3) The indirect hyperfine field is contributed by a 
macro scopically high number of QD nuclei. The strength 
of the indirect hyperfine interaction between nuclei in- 
creases for longer t c , and reaches a maximum for Ot c > 
1. This maximum is about huj^/fl. It is inversely pro- 
portional to nuclear polarization. Usually the field is 
less than the Knight field, hui n SQ. The indirect field 
plays an important role in the problem of electron spin 
dephasin g 17 ' 18 i 19 . It may also be important in the real- 
ization of dynamic nuclear self-polarization^, where elec- 
tron polarization and Knight field are equal to zero. 



FIG. 9: Dependence of nuclear spin relaxation time on nuclear 
spin polarization. Circles are for N = 251, squares for N = 
485, diamonds for N = 895, and triangles for N = 1365. 
All points coincide well with the straight line calculated from 
Eq. pi), with slope T 2 >y w 3. 



larization. T^(/3, N)(u) /N has approximately the same 
value for different N, and calculated points for different 
system sizes are grouped around a straight line. As fol- 
lows from Eq. (|23|) . the slope of this line gives the charac- 
teristic time, which in this case is T2{u>) ~ 60 or T27 ~ 3. 
These results are also in good agreement with Eq. (f!?3"|k 
they demonstrate the universal character of the predicted 
dependence of the dipole-dipole relaxation time for the 
nuclear spin potential £ on pn(/3) and system size, N. 



IV. SUMMARY 

(1) We demonstrated that for the short correlation time 
limit the rate of nuclear relaxation on electrons is propor- 
tional to the correlation time, i. e., T\ e cx T~ l . In the op- 
posite limit of long correlation time (£It c 3> I) we found 



(4) In the regime of long correlation time, the state of 
the quantum dot is characterized by three thermodynamic 
potentials: c, £ and %. These potentials directly affect 
the average electron spin, the average nuclear spin, and 
the nuclear spin temperature, respectively. The relax- 
ation of the nuclear spin potential £ by the dipole-dipole 
interaction affects the transfer of angular momentum to 
the crystal lattice. The diffusion of energy to the QD's 
environment is the main mechanism for the relaxation of 
the inverse nuclear spin temperature, X- The relaxation 
of the electron spin potential <; is connected with phonon 
scattering. 

(5) Our numerical simulation showed the ENSS's be- 
havior for the simple case of a pure hyperfine interaction 
between a resident electron and nuclei, and for the real 
case of an additional dipole-dipole interaction between 
nearest nuclei. They demonstrated a suppression of the 
£ relaxation that was caused by a decrease in the nuclear 
spin temperature. The dipole-dipole relaxation time was 
found to be proportional to the number of nuclei, and 
to the nuclear polarization squared. These results arc 
in good agreement with the analytical expression intro- 
duced in Refs.HHS 
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where (•■■)« represents the average over the initial spin 
direction. The time dependence of electron spin is given 
by Eq. ((6|). Eq. (|A.2|) is an odd function of the electron 
spin, and to linear approximation it does not depend on 
polarization. For a random initial electron spin distribu- 
tion (s a {0)sp(0)) = S aj p\\s\\/3 = 5 a ^/A, yielding 



1 

Tr 







-(1 + 

4 I 1 



s(t) ■ J s(t)cIt^J exp /- 
2 



— >dTr.) = 



(Or c ) : 



(A.3) 



APPENDIX: DERIVATION OF NUCLEAR SPIN 
PRECESSION EQUATIONS 

To derive Eq. JTJJI) we introduce Eq. (TO]) in Eq. © 



dt 



«o- 



(»t c ) 2 s - [ft x gg] ■ t c - (ft ■ sp) ■ fir 2 

1 + (fir c ) 2 

(A.l) 

For small nuclear polarization |(J)| <g; ||/||, and the vec- 
tor product in the right-hand side of Eq. (jA.lj) plays the 
main role in the average rate of nuclear polarization, 
i. e., ((fir c ) 2 \s x I n 



[ft X S ] X I„])/T c 



)/, where ((fts ) [ft x I„] r 2 )/ < 
-|cj n r c ||/|| 2 s . Here (■■•)/ repre- 
sents an average over the nuclear spin direction. Eq. (JTUJ) 
directly follows from this. 
It follows from (Eq. [12]) that 

(*fc). = £ </„" {(•(*) ■ In) Si s(r)dr 

- I n • /* s(r)dr) } exp{-t/r c }dt) a , (A.2) 



and 



[ft X I n ] T 



and finally 



i™±2 



1 + (Qt c ) 2 J l + (r>r c ) 2 



(A.4) 



/din 
\ dt 



2J„ 



(fir c ) 2 ) 



1 + (r!r c ) 2 



[ft X I n ]T c 



l + (^r c ) 2 J- (A - 5) 
The first term in the right-hand side of Eq. (|A.5[) con- 
tains the relaxation of the nuclear polarization compo- 
nents both parallel (J n ||) to and transverse (I n ±) to ft, 
see Eq. (|T2l . Its last term describes the nuclear spin 
precession in the indirect hyperfine field, see Eq. (|13l) . 
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||t/>(-R„)|| 2 and ui n depend on the type of nucleus and on 



its position in the crystal cells and QD. 

In this paper the field is given in energy units. The negative 

sign of the Knight field implies that the hyperfine energy 

has a minimum when I and so are parallel. 

We neglect the part of Eq. ((9} proportional to the square 

of the mean electron spin. 

This saturation value is determined by the fluctuation 
of the electron spin along SI. The average values of nu- 
clear spins, (I), and frequency, (SI), are parallel to each 
other, and the relaxation rate for the polarized nuclei, 
7|j > \\I\\/VN, decreases as (1 + (SlTc) 2 )^ 1 . For states 
with low polarization the fluctuation of the nuclear spin 
directions are more important than the fluctuation of the 
modulus. For these states, increasing the correlation time 
decreases relaxation rate by a factor of 2. 
A fluctuating magnetic polaron was investigated in the Ra- 
man spectros copy of diluted magnetic semiconductors in 

Refs. [alilUliJ 

As is well known from magnetism theory, the energy of 
magnetic anisotropy for a cubic spin lattice is proportional 
1° Is x + Is y + Is,z + 0(/|). It is important only for high 
enough nuclear polarization 



